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The origin of the charge ordering in organic compounds ^-(BEDT-TTF)2X {X = 
MM'(SCN)4, M =Tl,Rb,Co, M' =Cs,Zn) is studied using an extended Hubbard model. Cal- 
culating the charge susceptibility within random phase approximation (RPA) , we find that the 
(3 X 3) ~ (3 X 4) charge ordering observed at relatively high temperatures can be considered 
as a consequence of a cooperation between the Fermi surface nesting, controlled by the hop- 
ping integral in the c direction, and the electron-electron interactions, where the distant (next 
nearest neighbor) interactions that have not been taken into account in most of the previous 
studies play an important role. Mean field analysis at T = also supports the RPA results, 
and further shows that in the 3x3 charge ordered state, some portions of the Fermi surface 
remain ungapped and are nested with a nesting vector close to the modulation wave vector of 
the horizontal stripe ordering observed at low temperatures in X = MM'(SCN)4. We further 
study the possibility of superconductivity by taking into account the distant off-site repulsions 
and the band structure corresponding to X in which superconductivity is experimentally 
observed. We find that there is a close competition between d^iy-wave-like singlet pairing and 
P2!+2j/-wave-like triplet pairing due to a cooperation between the charge and the spin fluc- 
tuations. The present analysis provides a possible unified understanding of the experimental 
phase diagram of the ^-(BEDT-TTF)2X family, ranging from a charge ordered insulator to a 
superconductor. 
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1. Introduction 

6'-(BEDT-TTF)2A: is one of the most interesting fam- 
ilies of organic compounds, ranging from charge ordered 
insulators such as X =RbZn(SCN)4 to a superconduc- 
tor X =13.-'^'^ The anion X controls the angle between 
the BEDT-TTF molecules, which in turn determines the 
band structure. A further fascination for this series of 
compounds has arisen by a resent observation of a giant 
nonhnear transport in X =CsA/'(SCN)4 (M' =Co,Zn), 
which makes the material work as an organic thyristor.'^' ^ 
There, it has been pointed out that the coexistence of 
two kinds of short range charge ordering with modula- 
tion wave vectors (Q'a,9b,'Zc) = (fj^ij) ^"^^ (O^^i^) (i^^ 
units of the reciprocal lattice primitive vectors in a, h, c 
directions; a and c directions are shown in Fig.l) plays an 
important role in this phenomenon. The (|, fc, i) modu- 
lation corresponds to a 3 x 3 ordering with respect to the 
a-c unit cell shown in Fig.l, while a charge modulation 
with (0, fc, |) is often referred to as the horizontal stripe 
ordering (Fig. 2). At high temperatures, a diffuse X-ray 
spot is observed at (|, A:, \),^'^ while the (0, k, \) struc- 
ture develops as temperature is lowered, and the system 
becomes more insulating. When the electric field is ap- 
plied, the (0, fc, i) horizontal stripe ordering is degraded, 
resulting in a recovery of the metallic behavior and thus 
the nonlinear transport. Related to this is an observation 
of lattice modulation with a wave vector (|, fc, 0.29 ^ i) 
in X =CsCo(SCN)4 under pressure of lOkbar, which has 
been attributed to a pressure induced 2kp charge den- 



sity wave (CDW) because the modulation wave vector 
coincides with the nesting vector of the Fermi surface.^ 
In X =RbAf'(SCN)4 {M' =Co,Zn) also, two kinds 
of charge ordering are involved. At high temperatures 
in the metallic phase, diffuse X-ray spots are observed 
at (|,A:,i) and (|,A:, |), which suggests the presence 
of short range 3x4 charge ordering.^' ^ Some anoma- 
lies are also observed in the NMR experiments in the 
similar temperature range. ^^"'^^ At around 200K, the sys- 
tem undergoes a metal-insulator transition, accompa- 
nied by a structural phase transition into the so-called 
6d phase, in which the unit cell is doubled in the c 
direction.^' In the 6d phase. X-ray diffraction 
measurements have revealed that a long range horizontal 
stripe charge ordering with the modulation wave vector 
(0,0, i) takes place.i-^''^'^ 

Theoretically, the origin of these charge orderings has 
been an issue of great interest. Although some order- 
ings are accompanied by structural phase transition or 
lattice modulation such as the horizontal stripe order- 
ing in X =RbA/'(SCN)4 and the (|, A:, 0.29 - \) or- 
dering (2fcF CDW) in X =CsM'(SCN)4 under pres- 
sure, such a lattice modulation is not observed when the 
(0, k, i) (horizontal stripe) short range order develops in 
X =CsM'(SCN)4, nor when 3 x (3 4) ordering occurs 
mX = MM' (SCN)4 at ambient pressure. Therefore, it is 
reasonable to assume that the electronic degrees of free- 
dom takes the initiative in these charge orderings, and 
in some cases the lattice modulation occurs as a con- 
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Fig. 1. The lattice structure of the cation layer of (9-(BEDT- 
TTF)2X. Vp, Vc are the nearest neighbor interactions, while Va, 
Vq, V2c are the next nearest neighbor interactions. The a-c unit 
cell is the usual unit cell, while we can use the x-y unit cell to 
unfold the BZ. 
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Fig. 2. (Color online) Horizontal stripe, vertical stripe, diagonal 
stripe, and c-axis three-fold charge orderings are shown. The 
green (or dark) molecules are the charge rich ones. 



sequence of the electronic charge modulation. In under- 
standing the origin of the charge ordering in 0-(BEDT- 
TTF)2X from a purely electronic point of view, an impor- 
tant progress was made by Seo, who included the near- 
est neighbor off-site repulsions and in the model 
Hamiltonian (see Fig.l).^^ There have also been studies 
on a model that neglects V^,^^'^^ but an estimation by 
Mori^"^ shows that ~ K for the 6'-(BEDT-TTF) com- 
pounds. In ref.ll, the energy of various types of charge 



ordering patterns, such as horizontal, vertical, and diago- 
nal stripes (Fig. 2), was calculated within the mean field 
approximation. There, it was found, however, that the 
horizontal stripe charge ordering is not stabilized for the 
0-phase lattice structure, while it is stabilized in a cer- 
tain parameter regime in the 6^ lattice structure, which 
has a two-fold modulation in the c-axis direction. Thus, 
within this analysis, the lattice modulation seems to be 
the cause, not a consequence, of the electronic charge 
stripe formation. Exact diagonalization studies on these 
stripe orders have also been performed, but in the 
mean while, an analysis by Mori in the static limit re- 
vealed a possibility of non-stripe charge ordering shown 
in Fig. 2,^^ which has a three-fold periodicity in the c- 
axis direction. We will call this the "c-axis three-fold" 
charge ordering in the present paper. Later on, Kaneko 
and Ogata extended Seo's mean field study to take into 
account the possibility of this c-axis three-fold order- 
ing, where they found that this ordering does have a 
lower energy then the vertical, diagonal, and horizontal 
stripes when ^ Vc- In these studies, a possible relation 
between the c-axis three-fold ordering and the diffuse X- 
ray spots observed at (|,/c, i ^ j) was suggested. A 
variational Monte Carlo study by Watanabe and Ogata 
found a subtle competition between the c-axis three-fold 
ordered state and the diagonal or vertical stripe states, 
but quite recently, Hotta et al. used exact diagonalization 
for a spinless half-filled model, which can be considered 
as an effective model of the 0-type compounds in the 
large on-site U limit, and showed the presence of c-axis 
three-fold charge correlation for Vp ^ Vc-^^ 

In the present study, we will further confirm that the c- 
axis three- fold charge correlation is strong when Vp ^ Vc- 
Moreover, we point out that while the c-axis three-fold 
charge ordering is uniform in the a-axis direction, the 
diffuse X-ray spots are observed at (|,/e, | ~ j), which 
corresponds to orderings with a three-fold periodicity in 
the a-axis direction. We stress here that the common 
(and thus probably the most essential) feature in the 
high temperature charge ordering is the three-fold peri- 
odicity in the a-axis direction {qa = |), while the peri- 
odicity in the c-axis direction is sensitive to the anions 
and/or the pressure. Thus, we conclude that neither the 
low temperature horizontal stripe nor the high temper- 
ature 3 X (3 ~ 4) charge orderings can be understood 
within purely electronic models that consider only Vp and 
Vc (with Vp ~ Vc) as off-site repulsions. We argue that 
the consideration of more distant (next nearest neighbor) 
electron-electron interactions is crucial, and propose that 
the cooperation between the Fermi surface nesting and 
the electron-electron interactions including those distant 
ones is the origin of the charge orderings. 

Another aim of the present study is to investigate 
the origin of superconductivity in X =13.^^ Since this 
material sits in the vicinity of the (nearly) charge or- 
dered materials in the experimental phase diagram,^ it 
is interesting to investigate whether the charge fluctua- 
tions can give rise to the occurrence of superconductiv- 
ity. The mechanism of superconductivity in 0-type com- 
pounds has been studied theoretically in the past,^^'^^' 
but the distant interactions beyond Vp and Vc that are 
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Fig. 3. (Color online) The effective lattice structure using the x-y 
unit cell. We specify the hoppings and the interactions by c, p, 
a, g, 2c, which specify the position relative to the center (red) 
site. 



necessary to understand the charge ordering pattern for 
X = MM' (SCN)4 were not included there. Since it is 
natural to assume that the range of the off-site interac- 
tion does not change drastically with the change of the 
anion, we investigate the possibility of superconductivity 
in the vicinity of the charge ordering phase in a model 
that takes into account the distant interactions and the 
band structure of X 

2. Formulation 

2.1 Model 

In the lattice structure of 6'-(BEDT-TTF)2X shown in 
Fig.l, the direction of the molecules alternates along the 
a axis resulting in the rectangular a-c unit cell, but this 
alternation is irrelevant as far as the hopping integrals in 
the tight binding model are concerned. Therefore, we can 
take a unit cell (x-y) which is half the size of the usual 
unit cell, thereby obtaining the effective lattice structure 
shown in Fig. 3. The corresponding Brillouin zone (BZ) 
becomes unfolded, and the relation between the unfolded 
and the folded BZ is shown in Fig. 4. Here, the wave vec- 
tors are denoted in units of the reciprocal lattice primi- 
tive vectors, and the subscripts uf or / indicate whether 
the reciprocal vectors are those of the unfolded or the 
folded BZ. The points (|, A:, ^ ~ |)/ at which the diffuse 
X-ray spots are observed in the high temperature regime 
oi X ^ MM'(SCN)4 faU on the diagonal hue noted as 
ka — ^, which is a line that satisfies — ky = ^ in 
the unfolded BZ. It can be clearly seen that these po- 
sitions are different from the modulation wave vector of 
the c-axis three-fold ordering, (0, §)/ = (|, ^)uf- 

We consider an extended Hubbard model on the lattice 
shown in Fig. 3, where the Hamiltonian is given in the 
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Fig. 4. (Color online) Upper panel: The relation between the 
folded and the unfolded BZ is shown. Modulation wave vectors 
corresponding to various orderings are shown. The wave vectors 
are presented in units of the reciprocal lattice primitive vectors. 
The red line shows the Fermi surface for tc = —0.3. Lower panel: 
The BZ in the upper panel is squeezed in the kx = —ky direction, 
and than placed repeatedly. The trapezoid area surrounded by 
the red dashed line corresponds to the area at which the X-ray 
diffuse rods are observed in Fig. 4(a) of ref.7. 
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where c^^ creates an electron with spin a —1,1 at site 
i, and rii^ = cl^Cia. The band filling, the average num- 
ber of electrons per site (molecule), is fixed at n = 1.5 
in accord with the actual material. In Fig. 3, the letters 
p, c, a, • • • denote the relative positions with respect 
to the red site in the center, and we use these letters 
to specify the range of the hopping integrals and the 
off-site repulsive interactions. As for the kinetic energy 
terms, we consider hoppings tp and tc, where ip(— O.leV) 
is taken as the unit of energy. U is the on-site repul- 
sive interaction. As for the off-site repulsions, we con- 
sider, in addition to Vp and Vc (nearest neighbor inter- 
actions) that are usually taken into account, more dis- 
tant interactions Va, Vq, and V2c (next nearest neigh- 
bor interactions, see also Fig.l). Consideration of these 
interactions is supported by Mori's estimation showing 
that these distant interactions are not so small com- 
pared to Vp and Vc-^"^ Also, there is a recent trans- 
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port experiment which indicates that the repulsive in- 
teraction in 6'-(BEDT-TTF)2MZn(SCN)4 (M =Cs,Rb) 
is long ranged. In fact, the presence of such a distant 
interaction has also been pointed out theoretically for 
quasi-one-dimensional organic compounds (TMTSF)2X 
in the context of coexisting 2kp spin and 2kF charge den- 



sity waves 
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ity. 



29-34 



and spin triplet pairing superconductiv- 



So there is a possibility that long ranged nature 



of the electron-electron interaction may be common in 
these organic materials. Since the values of the electron- 
electron repulsions have not so far been strictly estimated 
from, e.g., first principles calculations, here we vary them 
as parameters and investigate whether the experimen- 
tally observed charge orderings can be understood within 
a realistic parameter regime, where we take Va,Vq, Vic to 
be smaller than ~ and 14 ~ to be smaller than 
~ VI2. 

2.2 Random Phase Approximation 

In this subsection, we describe the random phase ap- 
proximation (RPA) adopted in the present study. Within 
RPA,^^'^^'^^ the spin and charge susceptibilities, Xs and 
Xc, respectively, are given as 

xo(g) 



1 - Uxo{q) 

xo{q) 



(2) 



^-(«) l + {U + 2Viq))xo{qy 
Here xo is the bare susceptibility given by 

with £^ being the band dispersion given as 

^fe = 2tp[cos{kx) + cos(fcy)] -I- 2tc cos{kx + ky) (4) 

and f{ep) = l/(exp(gp — fJ.)/T) + 1) is the Fermi distri- 
bution function. When the nesting of the Fermi surface 
is good, Xoiq) is maximized at the nesting vector. V{q) 
is the Fourier transform of the off-site repulsions, given 
as 

V{q) = 2Vp[coa{q^) + cos{qy)] 

+ 2Vc cos{q^ + Qy) + 2Va cos{qx - Qy) 

+ 2Vq [cos(2(7j; +qy) + cos{qa, + 2qy)] 

+ 2V2cCOs{2q^ + 2qy) (5) 

To discuss superconductivity, the effective pairing in- 
teractions for the singlet and triplet channels due to spin 
and charge fluctuations are given as 

KZ%etiq) = u + v{q) + lu\M 



l{U + 2V{q)rxM 



.-{U + 2V{q)fxc{q). 



(6) 



(7) 



To obtain the superconducting transition temperature 
Tc, we solve the linearized gap equation within the weak- 
coupling theory, 

1 tanh{/3er.//2) 
AA(fc) = - fe') i^-fe^A(fe'). 

(8) 

The eigenfunction A of this eigenvalue equation is the 
gap function. The transition temperature Tc is deter- 
mined as the temperature where the largest eigenvalue 
A reaches unity. In the actual numerical calculations, we 
take up to A/' = 256 x 256 fc-point meshes. 

2.3 Mean Field Approximation 

In the mean field approximation, we approximate the 
interaction terms niaUja' by n^o- {nja' ) + {nia)njai . We as- 
sume an electron density modulation with a modulation 
vector (i, —\)uf (two-fold periodicity in the x direction 
and six- fold periodicity in the y direction), resulting in a 
unit cell that contains six sites (see the inset of Fig. 12 in 
section 3.2). This corresponds to the 3x3 ordering in the 
original lattice structure because (§,5)/ = {\,—\)uf- 
Here we neglect the spin ordering for simplicity, i.e., 
(nji") = (n^i), because (i) the ordering of the spins, if 
any, should take place in accord with the nesting vector 
of the Fermi surface, and therefore incommensurate and 
diSicult to deal with, and (ii) the main aim of the mean 
field analysis is just to reinforce the RPA results and also 
to obtain the Fermi surfaces in the charge ordered state. 
We do not compare the total energy between different or- 
dering states to pin down the ground state because (i) we 
do not consider the spin ordering properly and (ii) there 
could be a number of candidates for the ground state 
including incommensurate states. We obtain the band 
structure in the presence of the (^, — g)u/ ordering, and 
then calculate {rucr) at T = 0, which is substituted back 
into the mean field Hamiltonian, and the band structure 
is calculated again. These procedures are repeated until 
all {nia) are self-consistently determined. We define 

Anmax = "maxT " ('T-t)> (9) 

where rimaxt is the largest up(=down) spin electron den- 
sity among the six sites and (n|) is the average up spin 
electron density, 0.75 in the present case. In the actual 
numerical calculations, we take up to 200 x 200 unit cells. 

3. Cheirge fluctuations and ordering 

3.1 Random Phase Approximation 

First we consider the case when tc = and Va ~Vq — 
Vic = 0. In Fig. 5, we show the RPA charge susceptibility 
for f7 = 3 and Vc = 1.5, In the case ofVp = 1.5, there is a 
peak at (i, ^)m/, which corresponds to (0, |)/ in the orig- 
inal BZ, namely, the modulation wave vector of the c-axis 
three-fold ordering. This is due to a peak in —V{q) at 
q = (i, ^)uf as shown in Fig. 5. Namely, although xo('7) 
is maximized around (0, ±i)„/ and (±i,0)„/, the effect 
of —V{q) dominates in the denominator of the RPA for- 
mula of the charge susceptibility eq.(2), thereby minimiz- 
ing the denominator and maximizing Xc{q) at (^, ■^)uf- 
For Vp slightly smaller than Vc, which is more realistic for 
6'-(BEDT-TTF) compounds,^'^ the peak at (i, i)„/ stiff 
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Fig. 5. (Color online) Contour plots of the RPA result of the 
charge susceptibility Xcj the bare susceptibility xOi and the 
Fourier transform of the off-site interactions —V{q) are shown. 

= 0, U = 3, Vp = Vc = 1.5, Va = Vq = V2c = 0, T = 0.05, 
all in units of tp . The arrows represent the wave vector at which 
each function takes its maximum. 
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Fig. 6. (Color online) Plots similar to Fig. 5 except that Vp = 
1.3. Other parameters are the same: tc = 0, U = 3, Vc = 1.5, 
Va = Vq = V2c = 0. In Xc, the dashed arrow shows the wave 
vector at which a subdominant peak exists. The dashed arrow 
in —V{q) indicates that the wave vectors at which —V{q) takes 
large values extend toward (i, 0)uf compared to the case shown 
in Fig. 5 

remains (Fig. 6), but there also appears a subdominant 
peak around a wave vector somewhat close to (^,0)„/, 
i.e., the modulation vector of the diagonal stripe. The 
origin of this second peak becomes more clear if we fur- 
ther introduce a finite tc = —0.3, which roughly corre- 
sponds to the case of X =RbM'(SCN)4 (Af'=Zn, Co) 
or X =CsCo(SCN)4 under pressure.^ In this case, the 
nesting of the Fermi surface (see Fig. 4) becomes good 
around (|, j)f = (jj, — ^)u/, as can be seen from the 
bare susceptibility xo shown in Fig. 7. Now the peak near 
{^,0)uf has about the same height as that at (^, ■!)„/, 
which is due to a combination of Xo('7) that has large val- 
ues extending from the nesting position toward (i, 0)„/, 
and —V{q) in which large values extend to wave vectors 
closer to (^,0)„/ when Vp < Vc- 

Thus, in all of the cases studied above, the peak posi- 
tion of the charge susceptibility disagree with the exper- 
imental observations. We have looked into various other 
cases with Vc ^ Vp and Va = Vq = V2c = 0, where we 
always found that the charge susceptibility does not have 
a peak at (i, -i)„/ - (ii, -^)„/ = (|, 5 j)f, corre- 
sponding to 3 X (3 ~ 4) modulation, nor at {j, j)uf = 
(0, ^)f, which corresponds to the horizontal stripe mod- 
ulation. 

Thus the model that considers only U, Vp, and Vc as 
electron-electron repulsions fails to explain the experi- 
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Fig. 7. (Color online) Plots similar to Fig. 6 except that tc = 
—0.3. Other parameters are the same: U = 3, Vp = 1.3, Vc = 1.5, 
Va = Vq = V'ic = 0. The dashed arrow in Xo(9) indicates that 
the wave vectors at which xo (q) takes large values extend toward 
the end point of this arrow. 
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Fig. 8. (Color online) Plots similar to Fig. 7 except that Va = 0.3, 
Vq = 0.5, V2c = 0.7. Other parameters are the same: tc = —0.3, 
U = 3,Vp = 1.3, Vc = 1.5, T = 0.05. 

mentally observed charge orderings. Since we believe that 
the electronic degrees of freedom is responsible for these 
orderings, we now turn on the distant interactions Va, Vq, 
Vic- In Figs. 8 and 9(a), we show the charge susceptibility 
for [/ = 3 and J7 = 4 with tc = —0.3, maintaining the 
ratios of the off-site repulsions to \J roughly the same. In 
both cases we find a peak near the positions at which the 
diffuse X-ray spots are observed in the metallic phase of 
X = MM'(SCN)4. This is due to the fact that in these 
cases, — V^(g) is now broadly maximized in a region that 
includes the nesting vector position, as shown in Figs. 8 
and 9(a). Therefore, the denominator in eq.(2) is now 
minimized at a position close to the nesting vector, re- 
sulting in the charge susceptibility peak position that 
nearly coincides with the diffuse X-ray spot positions. 

In Fig. 9(b), we show the result for tc — —0.05 with 
other parameters being the same with (a). This choice 
of tc corresponds to X =CsM'(SCN)4 at ambient pres- 
sure.^ In this case, the peak position of the charge suscep- 
tibility moves toward the unfolded BZ edge, i.e., closer to 
(\^ ~\)ui — (|j as compared to the case with tc — 
—0.3. This is because the nesting vector of the Fermi sur- 
face is closer to (§,5)/ as can be seen in the figure of xo- 
This tendency is again consistent with the experimental 
fact that diffused X-ray spots at high temperatures are 
observed at = | in X =CsM'(SCN)4 at ambient pres- 
sure, while they are &i \ < qc <\\i\ X =CsM'(SCN)4 
under pressure and in X =RbM'(SCN)4. 

Here we have presented the charge susceptibility for 
certain choices of interaction values, but we have also 
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Fig. 9. (Color online) (a) Plots similar to Fig. 8 with tc = —0.3, 
U = 4, Vp = 1.8, Vc = 2.0, Va = 0.4, Vq = 0.7, V2c = 1.1, 
T = 0.25. (b) Plots similar to (a) except that tc = -0.05. V{q) 
is the same with that shown in (a). 



investigated cases with other sets of parameter values 
and found similar peak positions of the charge suscep- 
tibility within a certain parameter regime. In Fig. 10, 
we show the values of the distant interactions for which 
the charge susceptibility has a peak close to the the dif- 
fuse X-ray spot positions, where "close" means {qx, qy)uf 

1 

' 6- 



within the range of 0.45 < Qx < ^ and —j < Qy < 
As seen from this figure, we may say that good agree- 
ment with the experimental observation, as far as the 
peak position of the charge susceptibility is concerned, is 
obtained when Va + Vq is roughly constant for a given 
V2c- The peak height (the strength of the charge fluc- 
tuations) on the other hand becomes higher for larger 
Va and smaller Vq, which is because —V{q) around 
1 = ih^l)^f is larger in this case, so 

that the cooperation between the off-site repulsions and 
the Fermi surface nesting becomes more effective. 

Thus, in our present view, the 3 x (3 ~ 4) charge or- 
dering found in the high temperature metallic phase of 
X = MM' (SCN)4 occurs due to a cooperation between 
the Fermi surface nesting and the effect of the off-site 
interactions including not only Vc and Vp but also Va, 
Vq, V2c- Note that the term "charge ordering" is usually 
used as a real space ordering of charges that is caused by 
the repulsive electron-electron interactions (regardless of 
the shape of the Fermi surface), while a "charge density 
wave" is caused mainly by the nesting of the Fermi sur- 
face. Since the Fermi surface nesting and the electron 
repulsion cooperate in the present view, the charge or- 
dering in the metallic phase oi X — MM' (SCN)4 can 
be considered as a mixture of these two concepts. In this 
sense, the 3 x (3 ~ 4) charge ordering is conceptually dif- 
ferent from the c-axis three-fold charge ordering in that 
the latter has nothing to do with the Fermi surface nest- 
ing. We believe that chances of the diffuse X-ray spots 



Fig. 10. (Color online) Marked by circles (V2c = 1-1) or by trian- 
gles (V2c = 0.7) are the values of {Va,Vq) for which the charge 
susceptibility peak exists near the diffuse X-ray spot positions 
in the metallic phase (see text) for tc = —0.3, U = 4, Vp = 1.8, 
and Vc = 2.0. The red (or dark) [green (or light)] circles are the 
values of (Va, Vq) for which triplet [singlet] pairing superconduc- 
tivity dominates for tc = -1-0.2 (see sections. 1). For (Va, Vq) with 
uncolored (open) circles, charge ordering occurs above T = 0.2 
for tc = +0.2. 



accidentally coinciding with the nesting vector are very 
slim, and that the Fermi surface nesting must have some- 
thing to do with the 3 x (3 ~ 4) charge ordering, while on 
the other hand, we do believe that the electronic degrees 
of freedom, not the lattice, is the origin of this ordering. 

Now, to further show that the value of tc is impor- 
tant for the charge fluctuations and ordering, we show in 
Fig. 11 the charge ordering temperature Tco, defined as 
the temperature at which the RPA charge susceptibility 
diverges at a certain q, as a function of tc- Tco decreases 
with the increase of tc (from tc < to tc > 0), and 
sharply drops for tc > 0. In fact, the increase of tc induces 
several effects. First, the nesting of the Fermi surface 
is degraded (compare Fermi surfaces in Fig. 4 {tc < 0) 
and Fig. 15 {tc > 0)), secondly, the nesting vector de- 
viates from the region where —V{q) is maximized, and 
thirdly, the density of states near the Fermi level becomes 
smaller. ■^^ These factors cooperate to work destructively 
against charge ordering. 

This tc dependence of Tco is reminiscent of the phase 
diagram of the 6'-(BEDT-TTF) family obtained by Mori 
et al.,^ where charge ordering takes place for compounds 
having large negative tc such as X =RbM'(SCN)4, while 
it docs not for compounds having small or positive tc like 
X However, note that in Fig.ll, the charge order- 
ing takes place at Tco with a modulation wave vector 
close to the nesting vector of the Fermi surface (around 
(|, I ~ -j)/ for tc < 0), while the charge ordering tem- 
perature in the experimental phase diagram is those for 
the (0, fc, ^)/ ordering, i.e., the horizontal stripe state. 
In fact, although we have succeeded in understanding 
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11. (Color online) Tco is plotted as a function of tc for U ■ 
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O.leV ^ lOOOK is taken as the unit of temperature. The 
values of tc/tp for X = MM'(SCN)4 has been taken from table 
II in ref.6. As for X =I3, the value of tp is taken from ref.2 while 
that of tc is from ref.42. 



the high temperature 3 x (3 ~ 4) ordering, the horizon- 
tal (0, fc, ordering remains to be explained. Namely, 
there is no trace of enhanced charge fluctuations around 
{j, j)uf — (0, 5)/ in the charge susceptibility shown in 
Figs. 8 and 9. We will come back to this point in the latter 
part of the next section. 

3.2 Mean Field Approximation 

We now move on to the mean field analysis at T = 0. 
As mentioned in section 2.3, we assume 3x3 ordering 
shown in Fig. 12. The actual calculation is done for the 
ordering shown in the inset. We take values of the in- 
teractions to be slightly larger than those adopted in 
the RPA calculation, U = 5, Vp = 2.2, K = 2.5, 
{Va, Vq, V2c) = a(0.5, 0.8, 1.2), thereby roughly maintain- 
ing the ratios of the off-site repulsions to U when a — 1.^* 
a is a parameter introduced so as to investigate the effect 
of the distant interactions in a continuous manner. 

In Fig. 13 (a), we show An^ax defined by eq.(9) as a 
function of tc for a — 1. Here again, we find a similar 
tendency as found in the RPA calculation, namely, large 
negative tc induces the 3x3 charge ordering. In Fig. 13(b), 
we plot Arimax as functions of a for various values of tc- 
Here also, we find a tendency that the distant interac- 
tions induce the 3x3 charge ordering. For tc = —0.3, 
we find a finite Arimax solution even for a = 0, which is 
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Fig. 12. (Color online) 3x3 and 3x4 charge ordering patterns 
are shown. The dashed red rectangles represent the unit cell. In 
the upper inset of the 3x3 ordering, the corresponding pattern 
on the effective x-y square lattice, for which the actual mean 
field calculation was done, is displayed along with the unit cell 
that contains 6 sites. 




Fig. 13. (Color online) (a) Arimax obtained in the mean field 
calculation plotted as a function of tc- (b) An^nax for various 
choices of tc plotted as functions of the parameter a that controls 
the distant (next nearest neighbor) part of the off-site repulsions. 
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Fig. 14. (Color online) (a) Fermi surface of the 3x3 charge or- 
dered state obtained by the mean field approximation. U = 5.0, 
Vp = 2.2, Vc = 2.5, Va = 0.5, = 0.8, = 1-2, and tc = -0.3. 
The original folded and unfolded BZ are also shown. The dashed 
lines represent the boundaries of the new BZ. See text for the ex- 
planation of the vectors (1)~ (3). (b) The expected Fermi surface 
for the 3x4 charge ordered state (see text). 

probably because the nesting vector of the Fermi surface 
best matches the 3x3 ordering at this particular tc, but 
we believe that the 3x3 charge ordering state is not the 
ground state at least when a = considering our RPA 
results. (Note that we do not compare the total energy 
between different states as mentioned in section2.3. ) 

Now we turn to the Fermi surface in the charge or- 
dered state. In Fig. 14(a), we show the Fermi surface in 
the 3x3 charge ordered state for ic = —0.3 and a = 1. 
Note that the new BZ and the Fermi surface here are 
"unfolded" in the sense that the unit cell in the inset of 
Fig. 12 neglects the alternation of the molecules in the a- 
direction. (Therefore, the area of the new BZ is not 1/9 of 
the folded BZ but 1/6 of the unfolded BZ.) We find that 
some portions of the Fermi surface remain ungapped, so 
that the system should be metallic, which is consistent 
with the experiments, where the system remains metallic 
when diffuse X-ray spots are observed at (|, fc, i ~ i)/. 
What is even more interesting is that there is a nesting 
of the remaining Fermi surface with a nesting vector (3) , 
which is close to {\,\)uf = (0,^)/, i.e., the modula- 
tion vector of the horizontal stripe ordering. This "new 



nesting" originates from the fact that there are actually 
two different nesting vectors in the original system by 
symmetry, (+|,i - \)j and {-\, \ ^ j) f ■ The for- 
mer nesting is used to form the present 3x3 order- 
ing, while the latter nesting (vector (2) and (2)', where 
(2) — (2)' -I- (0, l)u/) remains in the ordered state, al- 
though modified to some extent from the original nesting 
vector because of the deformation of the Fermi surface. 
In Fig. 14(a), vector (1)=(-|-|, i)y is a reciprocal lattice 
vector of the 3x3 ordered state at which a finite Fourier 
component of the charge density exists, so that the re- 
maining portions of the original Fermi surface translated 
by vector (1) can interact with one another in the ordered 
state. Since (2) + (1) = (3) holds, the original nesting at 
vector (2) results in a "new nesting" of the Fermi sur- 
face at vector (3) in the ordered state, which is close to 
(0, i)/ especially when the c component of vector (2) 
(modified due to ordering) is relatively small. Since we 
know from the RPA results that the charge fluctuations 
are enhanced in the non-ordered state around vector (2) 
due to the electron-electron interactions, we may expect 
that the charge fluctuations are enhanced around the new 
nesting vector (3) in the 3x3 ordered state due to the 
electron-electron interactions, although this remains to 
be confirmed in a future study. 

Although our mean field analysis is restricted to the 
3x3 charge ordering, we can predict what the Fermi 
surface should look like in the case of 3 x 4 ordering from 
the analogy to the 3x3 case. In Fig. 14(b), we show the 
expected Fermi surface when the 3x4 ordering occurs, 
where the dashed pink lines are the portions that are 
expected to vanish due to the (+|, \)f nesting of the 
original Fermi surface. The remaining Fermi surface in 
the 3x4 ordered state is expected to be nested to some 
extent with the nesting vector close to (— |, \)f (vectors 
(2) and (2)', where (2) = (2)' + (0, 1)„/). In Fig.l4(b), 
(2)-f (1) = (3) holds, where vector (1) = (+|, \)f is a re- 
ciprocal lattice vector of the 3x4 ordered state at which 
a finite Fourier component of the charge density exists. 
The resulting "new nesting vector" (3) is precisely (0, ^)/ 
when vector (2) is (— |, \)f- This precise coincidence be- 
tween the new nesting vector and the horizontal stripe 
modulation vector comes from the fact that qc — j oi 
both the 3x4 modulation wave vector (1) and the nest- 
ing vector (2) is exactly half oi Qc — ^ of the horizontal 
stripe modulation vector. 

Now let us recall that within the RPA approach, we 
could not find any trace of (0, i)/ charge fluctuations. We 
propose here a possibility that the new ~ (0, |)/ Fermi 
surface nesting in the 3 x (3 ~ 4) charge ordered metallic 
state triggers, or at least plays some role in, the occur- 
rence of the horizontal (short or long ranged) charge or- 
dering in X — MM'(SCN)4 at low temperatures. In this 
view, the two charge orderings, the horizontal stripe and 
3 X (3 ~ 4), are not independent competing orderings, 
but the former is "based" on the latter, and therefore 
we may call the 3 x (3 4) ordering the "first stage or- 
dering" , and the horizontal stripe ordering the "second 
stage" . In this sense, the horizontal stripe should be more 
"fragile" as compared to 3 x (3 4) , which explains why 
it is the horizontal stripe ordering that is destroyed by 
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applying the electric field in X =CsAf '(SCN)4.'' Also, 
when the second stage ordering takes place, most of the 
remaining ungapped Fermi surface in the first stage or- 
dering should disappear, which also explains the (nearly) 
insulating behavior when the horizontal stripe ordering 
sets in. 

Let us now go back to Fig. 11 obtained within RPA. 
In the present terminology, T^^ is the temperature where 
the first stage charge ordering takes place. If the second 
stage charge ordering is indeed based on the first stage 
ordering, then we may consider that this tc dependence 
of Tco should also roughly represent the tc dependence of 
the temperature of the second stage ordering at (0, 5)/. 
In this sense, we may say that Fig. 11 is closely related to 
the experimental phase diagram by Mori et al.^ This view 
is also consistent with an experimental analysis under 
uniaxial pressure, which concludes that the value of tc is 
the key that dominates the occurrence of the horizontal 
charge ordering.**^ 

4. Possibility of Superconductivity: an RPA 
analysis 

Having found that the distant off-site interactions 
are necessary in order to understand the charge ordcr- 
ings and fluctuations of 6'-(BEDT-TTF) compounds with 

< 0, we now investigate the possibility of superconduc- 
tivity induced by such charge fluctuations. Since super- 
conductivity is observed in X =13, which has a positive 
^^42-44 concentrate on the case of tc > 0. In fact, 
the value tc/tp ~ 0.19 estimated by Mori for X =13,^^ 
falls close to a point at which Tco falls rapidly in Fig.ll, 
so that a charge fluctuation mediated superconductivity 
is possible. Here we mainly focus on two sets of off-site 
interaction values that give the correct charge suscepti- 
bility peak positions at tc ~ —0.3, i.e., U = 4, Vp — 1.8, 
Vc = 2 with Va = 0.4, Vq ^ 0.7, V2C = 1.1 or with 
Va = 0.6, Vq = 0.6, V2c = 1.1 (see Fig.lO). Using the 
RPA-fgap equation approach, we search for supercon- 
ductivity within the temperature range of T > 0.002tp 
under the condition that the denominator in the charge 
susceptibility formula (eq.(2)) is larger than 0.02. 

First, let us show the results for Va — 0.4, Vq — 0.7, 
V2C = 1.1. In this case, we find that the spin singlet pair- 
ing dominates over triplet pairing. The obtained spin- 
singlet gap function is shown in Fig. 15 (Agmgiot) along 
with the Fermi surface. The gap function changes sign 
as H 1 — along the Fermi surface with nodes inter- 
secting at positions on the and ky axes as in the d^y 
pairing. However, the entire gap function roughly has the 
form 

A(A;) oc cos(fc2: — fcj,) — 0.8 cos(fca; -I- fcj^) 

+ 0.3[cos(3fc2; -1- 2ky) + cos(2fca; + 3ky)] 

— 0.2[cos(2fcj; — ky) + cos{kx — 2fcj,)] 

- 0.2cos{4:k^+4:ky), (10) 
so that there is a deviation from the simple d^y gap form 








cos(fc2, 



cos{kx + ky) oc sin(fcj,) sin(fcj,). The addi- 



tional terms such as cos(3fca; -I- 2ky), cos(2fc2; — ky), and 
cos(4fc2, -|- 'iky) imply that the pairing occurs to some 



Fig. 15. (Color online) Contour plots of the spin and charge sus- 
ceptibilities, singlet and triplet pairing interactions, and the sin- 
glet {dxy-wave like) and triplet (pj;_|_2y-wave like) gap functions 
for U = 4.0, Vp = 1.8, Vc = 2.0, Va = 0.4, Vq = 0.7, V2c = 1.1, 
tc = 0.14, and T = 0.01. The gap functions are those that give 
the largest eigenvalue A for the singlet and triplet pairing chan- 
nels. The Fermi surface for tc = 0.14 (green solid line) is super- 
posed in the figures of the gap function. The dashed lines in the 
gap functions represent the nodal lines. Spin and charge fluc- 
tuations around wave vectors (1) ~ (3) play important roles in 
determining the pairing symmetry (see text). The yellow solid 
(light blue dashed) arrows indicate positive (negative) values of 
the pairing interactions. 



extent at distances {Ax, Ay) = ±(3, 2),±(2, — 1), and 
±(4,4) (in units of the lattice constants in the x and 
y directions), respectively, which can be considered as a 
consequence of the distant off-site repulsions that push 
away the electrons to far distances. 

Tc is plotted as a function of tc in Fig. 16(a), where it is 
enhanced in the vicinity of the charge ordering phase. For 
comparison, we also plot Tc in the absence of all the off- 
site repulsions, where the pairing is solely due to the spin 
fluctuations, and the charge fluctuations virtually do not 
contribute. Since Tc is higher in the presence of the off- 
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Fig. 16. (Color online) (a) Red solid line : Tc for the singlet dxy- 
wave pairing as a function of tc with U = 4.0, Vp = 1.8, Vc = 2.0, 
Va = 0.4, Vq = 0.7, and V2c = l-l- Green dashed line: Same as in 
the red solid line except all the off-site repulsions are turned off. 
(b) Tc for the triplet p2,+2j/-wave pairing as a function of tc with 
U = 4.0, Vp = 1.8, Vc = 2.0, Va = 0.6, Vq = 0.6, and V2c = 1.1. 
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Fig. 17. (Color online) (a) The eigenvalues of the gap equation 
for the spin singlet and triplet pairings, along with the ratio 
7 that measures the relative strength ratio between the spin 
and charge fluctuations (see text), are plotted as functions of /3, 
where the distant interactions are taken as (V^, Vc, Va, Vq, V2c) = 
/3(1.8,2.0,0.4,0.7, 1.1). U = 4.0, tc = 0.16, and T = 0.02 are 
taken. The yellow hatched area is the regime where charge fluctu- 
ations dominate over spin fluctuations, (b) Same with (a) except 
{Vp,Vc,Va,Vq,V2c) = /3(1.8,2.0,0.6,0.6, 1.1) and tc = 0.2. 
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site repulsions while the pairing symmetry is unchanged, 
we can say that the charge fluctuations cooperate with 
the spin fluctuations to enhance this superconductivity. 
This is also seen in Fig. 17(a), where we plot the max- 
imum eigenvalue of the gap equation A for the singlet 
and the triplet pairings for cases in which all the ofT-site 
repulsions {Vc,Vp,Va,Vq,V2c) are multiphed by a factor 
P iP = 1 corresponds to Vc = 2, Vp ^ 1.8, K = 0.4, 
Vq = 0.7, V2c = l-l)- Here we also plot a ratio 7 defined 
as 

_ [-{[/ + 2y(g)}xo(g)].nax , . 

^ [C/X0(9)]max ■ ^ ^ 

The denominator and the numerator in this ratio are 
terms that appear in the denominator of the spin and 
charge susceptibilities, respectively, and the subscript 
"max" denotes the maximum value within the entire BZ. 
Roughly speaking, the charge fluctuations dominate over 
spin fluctuations (the peak value of the charge suscepti- 
bility exceeds that of the spin susceptibility) when 7 ex- 
ceeds unity. We can see from Fig. 17(a) that in the charge 
fluctuation dominating regime, the spin singlet pairing 
(as well as the triplet pairing) is indeed enhanced as the 



Fig. 18. (Color online) Plots similar to Fig. 15 except that Va = 
0.6, Vq = 0.6, V2c = 1.1, and tc = 0.2. 

off-site repulsions are increased. 

To see this point in more detail, we plot the spin 
and charge susceptibilities, the singlet and triplet pair- 
ing interactions and the gap functions in Fig. 15. The 
spin susceptibility has a peak around wave vector (1), 
which bridges the portions of the Fermi surface with the 
same sign of the gap in the singlet d^y pairing, and there- 
fore work destructively against this pairing. (Note that 
A(fe)A(fc')yP"'^(fc-fc') < 0, fc, k' e Fermi surface, is fa- 
vorable for superconductivity with a gap function A(fc).) 
However, the spin susceptibility takes relatively large val- 
ues also around wave vectors (3)^(2), which bridges the 
portions of the Fermi surface with different gap signs 
in the dxy pairing. So the spin fluctuations have com- 
ponents which favor (around (1)) and unfavor (around 
(2) ^ (3)) the dxy pairing. On the other hand, the charge 
susceptibility strongly peaks near position (1), very close 
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to the peak position of the spin susceptibility, at which 
the spin fluctuations work destructively against the d^y 
pairing. The charge fluctuation around position (1) is 
so strong that it cancels out the spin fluctuation con- 
tribution (see eq.(6)) and even turns the singlet pairing 
interaction {V^^^gi^^ in Fig. 15) into a negative value at 
this position, thereby enhancing the dxy pairing. 

Now let us turn to the case oiVa = Vq = 0.6, V2c = l-l- 
Although the changes in the parameter values are not so 
large, and the peak position of the charge susceptibility 
for tc = —0.3 is very close to that in the case of Va = 0.4, 
Vq = 0.7, V2c = 1-1, we find that spin triplet pairing with 
a gap shown in Fig. 18 (Atripiet) dominates over singlet in 
this case. Tc of the triplet pairing is shown as a function 
of tc in Fig. 16(b), where superconductivity is again en- 
hanced in the vicinity of the charge ordering phase. The 
reason why spin-triplet pairing is enhanced is because the 
charge fluctuations arise around positions (1) ^ (2) (note 
that {^,ky) and {—^,ky) arc identical points) shown in 
Fig. 18, thus giving rise to attractive triplet pairing inter- 
action around position (2), which favors the gap Atripiet 
because vector (2) bridges the portions of the Fermi sur- 
face that have the same gap sign. An analysis similar to 
the one done for the singlet pairing shows that the gap 
roughly has the form 

A(fe) oc sm{kx + '2ky) — sin(2fcx + ky) 

- 0.3 sin(fc^ - ky) + 0.3[sin(2A;^) - sin(2fcy)] 

- 0.3[sin(4fcx) - sin(4fcy)] 

+ 0.3[sin(fc3, - 2ky) - sm{ky - 2ka:)]. (12) 

The dominating part has the form oc sm{kx + 2ky) — 
sm{2kx + ky), so we will call this pairing Pa;+23/-wave pair- 
ing. Although the pairing occurs mainly between sites 
separated by ±(1,2) and ±(2,1), here again, the pairs 
are also formed, to some extent, at far distances. The 
difference from the case of = 0.4, Vq = 0.7 is that the 
charge fluctuations tend to spread out toward position 
(2) in the case of Va = Vq = 0.6, so that the absolute 
value of the attractive triplet pairing interaction around 
position (2) is larger (~ 50 for Va ~ Vq ~ 0.6 and ~ 20 
for Va = 0.4, Vq = 0.7 In the case of K = 0.4, 
Vq = 0.7, the charge fluctuations are more localized 
around position (1), at which the charge fluctuations 
do not contribute to Pa;+2j/-wave pairing because vector 
(1) bridges the portion of the Fermi surface at which 
the nodes of the ]3j;+2y-wave gap exist. This difference 
of the charge susceptibility between the two cases origi- 
nates from the structure of V{q) discussed in section 3.1 
(note that the value of tc docs not affect V{q)), namely, 
-V{q) is larger around q = (|,-g)„/ ~ (ii'-^)"/' 
i.e., position (2), for smaller Vq and larger Va, so that 
the charge fluctuations spread out toward position (2) 
even when tc > 0, although the peak itself is at posi- 
tion (1) for tc > because the nesting vector is now 
close to that position. As a consequence of the differ- 
ence in the charge fluctuations, a difference arises also in 
the spin singlet pairing gap function. Namely, the strong 
charge fluctuations around position (2) in the case of 
Va = Vq = 0.6 results in attractive singlet pairing in- 



teraction not only near position (1) as in Va = 0.4, 
Vq = 0.7, but also near position (2), which works de- 
structively against dxy pairing, resulting in a vanishing 
gap around k (end point of vector (2) in Agingiet 

of Fig. 18) in the singlet channel, thereby giving way to 
spin-triplet pairing. In fact, we can see from Fig. 17 that 
triplet pairing dominates over singlet in the charge fluc- 
tuation dominating regime when {Va,Vq) = /3(0.6,0.6) 
(Fig. 17(b)), while singlet pairing continues to dominate 
deep within the charge fluctuation dominating regime 
when (K.Fg) =/3(0.4,0.7) (Fig.l7(a)). 

So the bottom line here for the superconductivity in 
the case of tc > is that spin triplet pairing tends to 
dominate over spin singlet pairing when the off-site in- 
teraction values are those that give strong charge fluc- 
tuations (namely high Tco) for < near the po- 
sition where the diffuse X-ray spots are observed in 
X = MM' (SCN)4. This tendency is more systematically 
seen in Fig. 10, where we classify the values of {Va,Vq) at 
which triplet (singlet) pairing dominates for tc = ±0.2 
by red (green) circles. Considering the ambiguity in the 
values of the on-site and off-site repulsions, we cannot 
predict at the present stage whether singlet or triplet 
pairings dominate in X =1^ . The singlet-triplet competi- 
tion is subtle, but in any case, the presence of the charge 
fluctuations enhances superconductivity. 

5. Comparison with (TMTSF)2X 

The coexistence of spin and charge fluctuations near 
the nesting vector "Q2fcjr" of the Fermi surface (al- 
though the nesting is not so good for tc > 0) is remi- 
niscent of a quasi-one-dimensional organic superconduc- 
tor (TMTSF)2PF6, in which we have proposed a pairing 
mechanism due to the coexistence of 2kp spin and 2kp 
charge fluctuations. ^^"^^ In this material, coexistence of 
2kF SDW and 2kF CDW has been experimentally ob- 
served in the insulating phase sitting next to the super- 
conducting phase in the pressure-temperature phase dia- 
gram.*^' This coexistence cannot be understood within 
a model that considers only the on-site U and the nearest 
neighbor repulsion V, and the effect of the second near- 
est neighbor repulsion V' is crucial,^^"^^'^^'^^'^^ sharing 
similarity with the present case. In fact, the mechanism 
of 2A;f(= f X 2 = f ) CDW in TMTSF is very similar 
to the mechanism of 3 x (3 ~ 4) ordering proposed in 
the present study. Namely, if wc consider only the near- 
est {V) and the next nearest neighbor (V) interactions 
within the chains, the Fourier transform of the off-site re- 
pulsions is given asy(q) = 2y cos((jj;) + 2y cos(2(72;)- On 
the other hand, the bare susceptibility Xo{q) has a peak 
at = ^ due to a good nesting of the Fermi surface of 
the quarter-filled quasi-lD band. Since both —V{q) and 
Xo{q) take large positive values at Qx = ^, 2V{q)xo{Q) 
in the denominator of eq.(2) takes a large negative value 
there, resulting in a peak in the charge susceptibility at 
2kF- Therefore, within this scenario, the 2kF CDW in 
(TMTSF)2PF6 occurs precisely due to a cooperation be- 
tween the effect of the off-site (next nearest neighbor) 
repulsions and the Fermi surface nesting. 

In the case of (TMTSF)2PF6, due to the good nest- 
ing of the Fermi surface and the one dimensionality, the 
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2kF peak positions as well as the overall q dependence 
of the spin and charge susceptibilities nearly coincide, so 
that the two fluctuations work destructively in the singlet 
pairing (cq.(6)) while constructively in the triplet pair- 
ing (eq.(7)). Also, due to the quasi-one-dimensionality of 
the system, the Fermi surface is disconnected in a man- 
ner that spin singlet d- wave-like pairing and spin-triplet 
/-wave-like pairing have the same number of nodes on 
the Fermi surface.^" Due to these factors favoring spin 
triplet pairing, charge fluctuations with similar strength 
with the spin fluctuation result in a subtle competition 
between /-wave and d-wave pairings in a realistic param- 
eter regime.^^'^^'^"' 

In the present case of the 6'-(BEDT-TTF) compounds, 
there are several differences as compared to TMTSF: 
(i) there is a large component in the spin fluctuations 
that work destructively against singlet pairing, (ii) the 
charge fluctuations cancel out this spin fluctuation com- 
ponent in the singlet pairing interaction to enhance the 
pairing, (iii) The singlet pairing has smaller number of 
nodes on the Fermi surface than in the triplet pairing 
because there is no disconnectivity in the Fermi surface, 
(iv) since there are 6 nearest neighbor {Vp, Vc) and 8 
next nearest neighbor (T4, Vq, V2c) interactions in the 6 
compounds as compared to 4 nearest (including the in- 
terchain direction) and 2 next nearest neighbor interac- 
tions in (TMTSF)2X, V{q) tends to be large, so that the 
charge fluctuations are Hkely to strongly dominate over 
spin fluctuations in a realistic parameter regime. In the 
charge fluctuation dominating regime, the charge fluctu- 
ations enhance both singlet and triplet pairings, but the 
enhancement in the latter is stronger. Factors (i) and (iv) 
favor triplet pairing, while (ii) and (iii) singlet pairing. 
Due to these competing factors, there is a close competi- 
tion between singlet and triplet pairing once again as in 
(TMTSF)2X. 

6. Conclusion 

In the present study, we have investigated the origin 
of the charge orderings in 6'-(BEDT-TTF)2AfAf'(SCN)4. 
We have shown that neither the low temperature hori- 
zontal stripe nor the high temperature 3 x (3 ^ 4) charge 
orderings can be understood within a model that takes 
into account only U, Vp, and Vc as electron-electron re- 
pulsions. We have shown that the 3 x (3 ~ 4) charge 
ordering can be considered as a consequence of the co- 
operation between the effect of the off-site repulsions in- 
cluding the distant interactions Va, Vq, and V2c, and the 
nesting of the Fermi surface. Moreover, we have proposed 
a possibility that the horizontal stripe charge ordering is 
triggered by the new (0, i)/ nesting in the 3 x (3 ~ 4) 
charge-ordered state, where some portions of the Fermi 
surface remain ungapped. From this viewpoint, we can 
understand why the horizontal stripe charge ordering 
temperature decreases with increasing tc (from tc < to 
tc > 0) in 6'-(BEDT-TTF)2X. Namely, as U increases, the 
nesting of the Fermi surface is degraded, the maximum 
positions of xo{q) and — V{q) deviate, and the density of 
states becomes small, all working destructively against 
3 X (3 ~ 4) ordering. Since 3 x (3 ~ 4) "first stage" or- 
dering gives base for the horizontal stripe "second stage" 



Author Name 

ordering in our view, the temperature of the latter order- 
ing should also become low as tc increases. 

In the present study, we have obtained the Fermi sur- 
face of the 3x3 charge ordered state by assuming that 
this first stage ordering is truly long ranged, while it is 
not a true long range order in the actual materials. Also, 
the role of the electron-electron interactions in the occur- 
rence of the second stage ordering as well as the second 
stage ordering temperature has not been addressed quan- 
titatively. A more detailed and quantitative understand- 
ing of this "successive charge ordering" view, as well as 
the investigation on its possible relevance to the nonlin- 
ear transport,"^' ^ serves as an interesting future study. 

Given that the distant interactions are important in 
6'-(BEDT-TTF)2X, we have further investigated the pos- 
sibility of superconductivity in the vicinity of the charge 
ordering phase for tc > 0. We have shown that there is 
a close competition between rfa;j/-wave-like singlet pair- 
ing and p2:+2y-wave-like triplet pairing, where the latter 
tends to dominate over the former when we adopt the 
interaction values that give strong charge fluctuations 
(high Tco) around (|, i j)f in the case of tc < 0. Re- 
gardless of whether triplet pairing dominates or not, the 
charge fluctuations are found to enhance superconduc- 
tivity. 

In total, we have provided a unified understanding of 
the overafl phase diagram of 6'-(BEDT-TTF)2X2 at least 
qualitatively within an electronic model that considers 
the band structure and the electron-electron interactions 
up to next nearest neighbors. 
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